LARGE TIME BEHAVIOR FOR VORTEX EVOLUTION IN 

THE HALF-PLANE 



D. Iftimie 
M.C. Lopes FilhoQ 

H.J. NUSSENZVEIG LOPESQ 



Abstract. In this article we study the long-time behavior of incompressible ideal flow in a half 
plane from the point of view of vortex scattering. Our main result is that certain asymptotic states 
for half-plane vortex dynamics decompose naturally into a nonlinear superposition of soliton-like 
states. Our approach is to combine techniques developed in the study of vortex confinement with 
weak convergence tools in order to study the asymptotic behavior of a self-similar rescaling of 
a solution of the incompressible 2D Euler equations on a half plane with compactly supported, 
nonnegative initial vorticity. 
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I. Introduction 

Let u = ou(t,x) be the vorticity associated to a solution of the incompressible two- 
dimensional Euler equations on the upper half-plane with an initial vorticity uoq which is 
bounded, compactly supported and nonnegative. We consider a rescaling W = W(t,x) = 
t 2 cj(t, tx), whose time-asymptotic behavior encodes information on the scattering of u into 

1 Research supported in part by CNPq grant #300.962/91-6 
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traveling wave solutions of the 2D Euler system on the half-plane. This choice of rescaling 
was also made in view of the fact that the horizontal velocity of the center of vorticity is 
bounded away from zero from below (see ||11|| ). The rescaling W is weakly compact as a 
time-dependent family of measures. The main purpose of this work is to present a structure 
theorem, stating that if the rescaling W is actually weakly convergent to a measure then 
this measure must be of the form ^ Tai5{x\ — aj) (g> S(x 2 ), with > 0, a discrete set of 
points on an interval of the form [0, M\ whose only possible accumulation point is X\ = 0, 
and where 5 denotes the one-dimensional Dirac measure centered at 0. 

Let us begin with a precise formulation of vortex dynamics on the half-plane. Let 
(x\,x 2 ) be the coordinates of a point x in the plane and denote the upper half-plane 
by H = {x 2 > 0}. The initial-boundary value problem for the incompressible 2D Euler 
equations in EI is given by: 



u t + u ■ Vco = 0, in (0, oo) x EI 

div u — 0, in [0, oo) x HI 

curl u = uj, in [0, oo) x H 

u 2 {t, x\,0) = 0, on [0, oo) x R 

uj(0, Xi, x 2 ) = ujq(xi, x 2 ) at {t = 0} x H, 



with u = (ui,u 2 
For bounded, 



the velocity and uj the vorticity of the flow. 

compactly supported initial vorticity, problem fll.lj ) is equivalent to the 
full plane problem with initial vorticity given by an odd extension of ujq to {x 2 < 0} (see 
FT] ] for details). Global well-posedness of the initial boundary value problem follows from 
this equivalence, using Yudovich's Theorem |37|. For compactly supported initial vorticity 



in (L 1 + BM + ) fl we have global existence of weak solutions by adapting Delort's 
Theorem to the half- plane case, see ||, [T7|, |S3|, |3"d] , but uniqueness is open. 

The present work is best understood within the context of research on vorticity confine- 
ment. Let uj = uj(t,x) be a (classical) solution of the full plane 2D Euler equations such 
that uj(0,x) is compactly supported. The problem of confinement of vorticity is to under- 
stand the spreading of the support of u(t, •) for large time. The main result in confinement 
of vorticity states that if the initial vorticity is nonnegative, with support contained in the 
ball -8(0; -Ro), then for any a > 1/4 there exists b > such that the support of vorticity 
at time t is contained in a ball of radius R(t) = (bt + ) a . This result is due to C. 

in lOl and 1134)1 . This area has seen 



Marchioro [18] for a = ^ and was improved to a > 



substantial recent activity, mostly in the direction of extending or generalizing Marchioro's 
original work, see |l|, [H| TT|, [16], [19, 



§ [H, ||, [23, |J 



Results on confinement of vorticity are rigorous actualizations of the rough idea that 
single signed 2D vorticity tends to rotate around, but not to spread out. This is false if the 
vorticity is not single signed, which can be seen by considering the behavior of vortex pairs, 
vorticity configurations that tend to translate to infinity with constant speed due to their 
self-induced velocity, see [JT3|] for a specific smooth example. Due to the traveling wave 
behavior of vortex pairs, vorticity scattering in two dimensions may become complicated, 
and interesting, when vorticity is allowed to change sign. Recently, the authors have 
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proved a new result on confinement of vorticity in this context. Let uj = uj(t,x) be a 
solution of incompressible 2D Euler in the full plane, with compactly supported but not 
necessarily single signed initial data u Q and let M — J u . For any a > 0, define the 
rescaling W(t, x) = t 2a uj(t,t a x). We have proved that if a > 1/2 then W(t, •) converges 
weakly to M5. This means confinement, in a weak sense, of the net vorticity in a region 
with roughly square-root in time growth in its diameter. This result will appear in a 



forthcoming article by the authors, [E2|. From the point of view of scattering, this result 
accounts for the behavior of the net vorticity, but says very little about the behavior of 
vortex pairs, because these tend to be weakly self-canceling when looked at from a large 
spatial scale. It one wants to study vortex scattering, the relevant information is the large- 
time behavior of \W(t, -)|, mainly in the case a = 1. The present article is directed precisely 
at this problem, with the simplifying assumption that the vorticity be odd with respect to 
a straight line, single-signed on each side of the symmetry line. Another way of expressing 
this is to say that in this article will study the scattering of co-axial, unidirectional vortex 
pairs. 

Let uj = uj(t,x) be the solution of the half-plane problem (|1 . 1| ) defined for all time, 
associated to initial data uj , which we assume, for simplicity, to be smooth, compactly 



supported and nonnegative. A confinement result proved by Iftimie in UTTJ, together with 
what we will prove here implies that the support of uj(t, •) is contained in a rectangle of the 
form (ai — bit a , ct) x (0, a<i + b%t^\ with a; real constants, bi, c > and < a, (3 < 1. We 
wish to examine the asymptotic behavior of the vorticity on the linearly growing horizontal 
scale that is naturally associated with the motion of vortex pairs. The approach we use 
is inspired on work on the asymptotic behavior of solutions of systems of conservation 
laws due to G. Q. Chen and H. Frid, see |1. Let u(t, x) = t 2 uj(t,tx). The function 
uj has bounded L 1 norm and will be shown to have support in a rectangle of the form 
(— bit a ~ l ,c) x (0,b 2 t l3 ~ 1 ). Hence the family of measures { u(t, -)}t>o is weak-* precompact 
and any weak limit of subsequences of this family is of the form fi®5o, with fi a nonnegative 
measure supported on the interval [0, c]. We will refer to such a measure fi as an asymptotic 
velocity density. Our main result may be stated in the following way. 



Theorem 1.1. Suppose that the initial data ujq for problem is such that there exists 
a unique asymptotic velocity density \i, i.e., uj(t, •) — /x<E> 5o when t — >• 00. Then ji is the 
sum of an at most countable set of Diracs whose supports may only accumulate at zero. 

The proof involves writing the PDE for the evolution of u and using the a priori es- 
timates available and the structure of the nonlinearity in a way that is characteristic of 
weak convergence methods, see 0. We will briefly discuss the physical meaning of both 
the hypothesis that u>(t, •) converges weakly and the conclusion regarding the structure of 
/i. 

The study of the wavelike behavior of vortex pairs goes back to Pocklington in p0[ , 
with more recent interest going back to work of Norbury, Deem and Zabusky and Pier- 
rehumbert, see ||, |26|, p9| . The existence (and abundance) of steady vortex pairs, which 
are traveling wave solution of the 2D incompressible Euler equations, i.e. vorticity shapes 
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which propagate with constant speed without deforming, has been established in the liter- 
ature in several ways, see p, [14], |26[. Steady vortex pairs have been object of an extensive 
literature, from asymptotic studies, see J36j and numerical studies, see [^TJ and even exper- 



imental work, see ||. Although some analytical results (see |[25|| ) and numerical evidence, 
pTfl , point to the orbital stability of steady vortex pairs under appropriate conditions, this 
stability is an interesting, largely open problem, see p2~ . 



Compactly supported vortex pairs interact in a way such that the intensity of the inter- 
action decays with the inverse of the square of the distance between them. Hence, vortex 
pairs moving with different speeds tend to behave like individual particles, decoupling af- 
ter a large time. This is what makes the study of vortex scattering interesting in this 
context. Let us illustrate the point of view we want to pursue with the example of the 
Korteweg-deVries equation. Nonlinear scattering for the KdV is well-understood, as solu- 
tions of KdV with smooth, compactly supported initial data are expected to resolve into 
a scattering state composed of an iV-soliton plus a slowly decaying dispersive tail. This 
fact was first formulated as a conjecture by P. Lax in |15| and broadly explored through 
the method of inverse scattering since then. The conclusion of Theorem [O] may be re- 
garded as a weak, or averaged form of Lax's conjecture for vortex pair dynamics. Note 
that steady vortex pairs correspond to classical solitons in this analogy, but no existence 
for the multibump solutions that would be associated to the classical N-solitons has been 
rigorously established. 

Let us call shape space the space of smooth compactly supported vorticity configura- 
tions, identifying configurations which are related through horizontal translations. Steady 
vortex pairs correspond to stationary shapes with respect to Euler dynamics. There are so- 
lutions of the two-dimensional incompressible Euler equations that describe periodic loops 
in shape space. Two examples of this behavior are: 1) a pair of like-signed point vortices 
on a half plane, which orbit one another periodically as they translate horizontally, called 
leapfrogging pairs, and 2) Deem and Zabusky's translational V-states, which are vortex 
patches with discrete symmetry, see H . From the point of view of scattering such solutions 
represent another kind of asymptotic state or, in other words, another kind of particle. Fur- 
thermore, one may well imagine solutions with quasiperiodic or chaotic behavior in shape 
space. Although there is no example of either case in the literature, the passive tracer 
dynamics of the leapfrogging pair is known to be chaotic, see [pgfl . Possible chaotic shapes 



represent an interesting illustration of Theorem |1.1| , as both the hypothesis of weak con- 



vergence and the conclusion are clearly related to the ergodicity of shape dynamics and 
the self-averaging of the velocity of the center of vorticity of such generalized vortex pairs. 
Finally, we must mention the work of Overman and Zabusky f2? |, where they do numerical 
experiments on the short term scattering of pairs of translational ^-states, the first (and 
only) study to date on the interaction of coaxial vortex pairs, which is the main point of 
the present work. 

The remainder of this article is divided into two large sections. The first one contains 
a discussion of confinement of vorticity in half plane flow, including two new results, hor- 
izontal confinement on the left and bounds on the distance to the boundary. The second 
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one contains the discussion leading to our main result, together with its proof. After this 
second section we include a brief section with conclusions and an Appendix containing a 
simple illustration of half-plane vortex dynamics. 



2. Confinement of vorticity 

2.1. Preliminary results. Let us begin by fixing basic notation. We denote by EI the 
horizontal half-plane given by HI = {x G M 2 ; x 2 > 0}. Reflection with respect to X2 = will 
be denoted by x = (xi,x 2 ) ^x = (xi, — x 2 ). If z = (zi,z 2 ) then its perpendicular vector 
is z L = (—z 2 ,Zi). We use L£(H) to denote the Lebesgue space of p-ih power integrable 
functions, p > 1, with compact support in HI. The dual of LP is LP , with the conjugate 
exponent given by p' = p/(p — 1)- The space of bounded Radon measures is denoted by 
BM and the Dirac delta at the origin is 5q. 

Consider the initial-boundary value problem for the incompressible 2D Euler equations 
in the half-plane fll.lp with initial vorticity Uq. If u is bounded then ( |1.1| ) is globally 
well-posed since it is equivalent, through the method of images, to an initial- value problem 
in the full-plane, with bounded, compactly supported initial vorticity (shown to be well- 
posed by Yudovich in [37fl). The method of images consists of the observation that the Euler 



equations are covariant with respect to mirror-symmetry. Thus an initial vorticity which 
is odd with respect to reflection about the horizontal axis will remain so, and give rise to 
flow under which the half-plane is invariant. Conversely, the odd extension, with respect 
to x 2 = 0, of vorticity in half-plane flow gives rise to full-plane flow. This observation 
is especially useful in order to deduce the Biot-Savart law for half-plane flow, to recover 
velocity from vorticity. 

Let us fix an initial vorticity uj q . We will assume throughout this article that ujq is a given 
nonnegative function in L p c {&) for some p > 2. If uj e Z^(HT), 2 < p < oo, then there exists 
a weak solution u, u> of (|1.1|) associated with this initial vorticity (see [PH ). Furthermore, 
u(t, •) > 0, t > 0, and the L 1 and L p -norms of u>(t, •) are bounded by the L 1 and L p -norms, 
respectively, of the initial vorticity. Using the method of images we can write the velocity 
u in terms of vorticity uj as: 



2-1) «(*,*)=/ n-,_ .12 ~ o I _-|2 "(t,V)dV. 



(x - y) L _ (x - y) 
-2tt\x — y\ 2 2tt\x — y 

We denote the kernel appearing the integral above by: 

{x-y) L (x-y) 



[2.2) K = K(x,y) 



2tt\x — y\ 2 27T |x — 2/| 2 ' 
whose components are given explicitly by: 

/ 0Q ^ v( \ y2[y 2 -x 2 2 + (xi -yi) 2 } ( , 2(x 1 -y 1 )x 2 y 2 

(2.3) K 1 (x,y) = : - 2 and K 2 [x,y) = — r=. 

7i\x — y\ z \x — y\ z 7r\x — y\\x — yy 
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It is easy to see that 
(2.4) \K(x,y)\< ' 



n\x — y\ 

from which we can deduce the fact that, if p > 2, then an L l fl L p -vorticity u gives rise to 
an L°°-velocity u with the estimate: 

ii n ^ n\\ up'/ 2 ii ll 1 -?'/ 2 

|F|U°°(H) < L,\\u\\ LPm \\uj\\ L1{m) . 

One interesting fact is that this estimate can be localized. From a technical point of view, 
this fact is the heart of the arguments presented in this work. To be more precise we will 
be using the following lemma, whose proof can be found in [11] : 

Lemma 2.1. Let a e (0, 2), S'cl 2 and let h : S —>■ R ; h > 0, be a function belonging to 
L 1 (S') f)L p (S), p > Then there exists a constant D = D a ^ p > such that the following 
estimate holds: 

I i7% dy ^ ^ll^l"(s)^llt(T) 72 ' v - e R2 - 
j s f — y\ 

Remark 2.1. To illustrate our use of this Lemma note that, when the function h is the 
vorticity and a = 1 this estimate implies that the portion of velocity due to vorticity in a 
region S will be small if the mass of vorticity in that region is small. 



2.2. One-sided horizontal confinement of vorticity. Iftimie showed in [11] that the 
horizontal component of the center of mass in half-plane flow travels with speed bounded 
below by a positive constant. This excludes any possible sublinear-in-time horizontal con- 
finement, at least in the direction xi > 0. On the other hand, half-plane flows with 
nonnegative vorticity have a tendency to move to the right, resisting left "back flow". The 
purpose of this section is to make this statement more precise. Our main result in this 
section is: 

Theorem 2.1. Let u 6 L%($H), p > 2, ujq > 0. Let u and uj be solutions of ( \L.lJ with 
initial vorticity uj$. Then there exists a positive constant D depending solely on the initial 
vorticity such that 

suppu(t, •) C {x e H ; xi > -D(tlogt)^} 

for all t > 2. 



Before we give the proof of Theorem |2.1| we need a technical lemma, in which we obtain 
an estimate on the mass of vorticity in the "back flow" region; we see that it is exponentially 
small. 

Lemma 2.2. Given k e N, there exist positive constants D\ and D2, depending only on 
the initial vorticity and on k, such that 

u{t,y) dy < — 

yi<—r r 



(i 



provided that r > D 2 (tlogt)z and t > 2. 



Proof. Consider the auxiliary function 77 = 77(5) 
ative, increasing and 



l+e s 



It is easy to see that 77 is nonneg- 



(2.5) 
Set 



\rj"(s)\ <rj(s) 



fr{t)= I rj(- Xl ^ r )u(t,x)dx, 



where A > will be chosen later. As 77 is nonnegative and increasing we clearly have: 



(2.6) 



fr(t) > 



xi<—r 



T]( ) u(t, x) dx > 77(0) 

Ar 



u(t, x) dx, 



x±<—r 



where we have used that for x\ < —r we have that 



Ar 



> 0. Therefore it suffices for our 



purposes to estimate f r (t). 

We will deduce a differential inequality for f r from which we estimate f r . To this end 
we differentiate in time to find: 

f' r (t) = -^- J ri'(- Xl ^ T )u 1 (t,x)uj(t,x) dx, 

where we have used the vorticity equation (|1 . 1| ) and integration by parts to throw deriva- 
tives onto 7], 



1 



2?rAr 



,/ x l +r. 



x 2 + V2 x 2 - y 2 



l \x — y\ 2 



\x — y\ 1 



u(t, x) u{t, y) dxdy, 



using the Biot-Savart law (271) 

1 

< 



2vrAr 



// x l +r.x 2 -y 2 ( v * v , , 
77 ' Ar - ' U _ ,,12 ^v*' x > V) dx d V> 



\x - y\ 



as 77', x 2 and 1/2 are positive. Finally, we symmetrize the kernel above by making the change 
of variables x y to obtain: 



1 



4n\r 



V ( T ) -V{- 

v Ar ' v 



Vi + r 
Ar 



x 2 - y 2 

N - y\ 2 



uj{t, x) uj{t, y) dx dy 



< — 1 1 — \v"(^x, y )\Tz — 7W x) cu(t, y) dx dy, 



AnXr 



Ar 



\x - y\ 



by the mean value theorem, with 9 x>y some point between and — 

Next we use Q2.5|) and the fact that 77 is nonnegative and increasing to deduce that 

|»m,v)l < \v(0 x , y )\<v(-^)+v(- y -^)- 



Ar 
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Since \x± — yi\ \x% — 2/2 1 < \x — y\ 2 we finally obtain the differential inequality: 



/;(*) < 



47rA 2 r 2 



7 ?(~ gl \ +7 ' ) +^(~ y \ +r )1 u(t,x)cj(t,y) dxdy = ^° 9 ^ / r (t), 



Ar ; ' ' v Ar ; J v ' ' v ' ^ ' 27rA 2 r 2 
where we have used that the L 1 -norm of u;(i, •) is constant in time. Integration now yields 

/,(»)< «0)«p(*^). 
Clearly we may assume, without loss of generality, that suppt^o C {xi > 0}. Then 

v( ^— )wo(x) dx < rj(—r) \\uq\\ia < exp(--) \\u Q \\ L i. 

Hence, we infer that 

/ p (*)<| Wxl «p(*^-i). 

In view of (|2.6| ), to finish the proof it is now sufficient to choose A such that 

CXP ('L A. 



— ) < — = exp(— klogr). 



The choice 



A 



2k logr 

is convenient provided that the following inequality holds 



(2.7) 



r 2 > ^ 2fc||^ || L i 



logr 



7T 



Notice that the function r i— > r 2 /logr is nondecreasing if r > e. Hence, choosing D2 
sufficiently large, it is easy to ensure ( |2.7|) if r > D 2 (tlogt)2 and t > 2. This completes 
the proof. □ 



Next we use Lemma 2.2 to estimate the horizontal velocity. 



Proposition 2.1. Under the hypothesis of Theorem \2. 1\ there exist positive constants D% 
and D4 such that 



\ui(t,x)\ < 



El 

\ x i\ 



for all t > 2 and x G EI such that X\ < — D^tlogt) 2 . 
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Proof. We will estimate directly i*i(t, x). From the Biot-Savart law ( |2.1| ) and the decay 
estimate ( |2.4| ) it follows that 

\v>i(t,x)\< [ -j— — :u(t,y)dy 



ir\x — y\ 



< ; H rc^(t,2/) d2/ + 

•yi<xi/2 n \ X V\ Jyi>xi/2 7l \ x y\ 



1 



|| P '/2 



7T 



I LP 



!/i<xi/2 7 ^Fl 



u(t,y)dy 
"ll^olUi, 



using Lemma [2.1| with a = 1. We have also used that both the L l and the L p -norms of 
u(t, •) are bounded by their initial values. 
Let 



k 



P 



+ 1, 



where [a] denotes the largest integer smaller than a. Choose D 2 as in Lemma |2.2| and let 
x satisfy x\ < — D^(t\ogt)^ with D4 = 2D 2 . The conclusion then follows from Lemma |2"T2 
with D 3 computed accordingly. □ 

We will finish this section with the proof of the horizontal confinement to the left. 



Proof of Theorem 2.1 . Let D% and D± be as in Proposition 271. If need be increase the val- 
ues of D 3 and _D 4 to show that any trajectory which reaches the region {xi < —D±{t log t) a } 
does not have enough horizontal velocity to go past the line x\ = — 2D 4 (tlogt)2. This 
proves that every trajectory lies in the region {xi > — 2D^{t logi) 2 } (with D± depending 
on the initial position of the trajectory); in particular, the support of the evolved vorticity 
stays in that region. □ 



2.3. Proximity to the boundary. For smooth flows it is easy to see that, if the support 
of the initial vorticity lies in the interior of the half-plane, then the support never reaches 
the boundary due to uniqueness of solutions of ODEs. Even if the flow is not smooth 
this remains accurate when understood in the context of the R. DiPerna and P.-L. Lions' 
theory of ODEs in Sobolev spaces J7|. However, we are left with no information on how far 
a particle path must remain from the boundary. The purpose of this section is to examine 
this issue. We will show that particle paths must stay at least as far as e~ CieC2± away 
from the boundary, for some positive constants Ci, C^. This result is interesting from the 
point of view of confinement, but it will not be used for the development of large time 
asymptotics. A similar concern was addressed by N. Depauw in his work on vortex patches 
in a bounded domain, see . 

We will use throughout this section the notation log to denote the map s 1— > log(s) = 

yl + log 2 (s). Note that | ^ logi | < \ for all t > 0. We denote a particle path by X = X(t), 
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so that 



jX=u(t,X). 

Theorem 2.2. Assume that ujq G L 1 D T/ien i/iere exists a positive constant D, 

depending only on the initial vorticity, such that every trajectory X verifies 



X 2 {t) > e 



-log(X 2 (0))e z 



for all times t > 0. 



Remark 2.2. If the initial vorticity is compactly supported in the half-plane then, as vortic- 
ity is transported by the flow, it follows that there exists a positive constant D\, depending 
only on the initial vorticity, such that 



support, •) C {x E M ; x 2 > e~ Die x *} 



for all times t > 0. 



Proof of Theorem \2.Q . We will prove that for all i6i, one has that 
(2.8) \u 2 (t,x)\<Dx 2 hgx 2 

for some constant D > depending only on the initial vorticity. Assuming that this is 
true, if X = (Xi,X 2 ) is a particle trajectory then: 



\X' 2 \ <DX 2 logX 2 , 



which in turn implies that 



— log log X 2 
at 



1 



logX 2 



X' 2 \og (X 2 ] 



< 



\X'\ 



< D. 



X 2 log Xo 



After integration, we infer that 

logbgX 2 < Dt + logbg(X 2 (0)), 

that is, 

k5gX 2 <l^g(X 2 (0))e D *. 

We therefore find that 

-logA 2 <bg(A 2 (0))e^, 

so that 

X 2 (t)>e-^ X2(0 ^ Dt , 



which proves Theorem 2.2. 
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Let us turn to the proof of ( |2.8| ). Fix x G H and start by noting that, by the Biot-Savart 
law Q) and (p|), 



2X2 

u 2 (t,x) = 



— — Vl V2 2 u(t,y)dy = ^ / Lx(x,y)uj(t,y)dy. 



7T J |S — y| 2 |x — 



7T 



We begin by observing that |xi — 2/i| < |x — y\ and, since x 2 , y 2 > 0, |y 2 | < |#2 + 2/2 1 < 
|5? — so that 



l-kifoy)! < 



a; — 2/| \x — 2/| 



It can be easily checked that 



x 2 + 2/2 > 



F2 - 2/2 1 + £2 



Using this estimate we find: 

,_ 1 . 1 n x . \xi -2/1I + \x 2 — 2/2 1 +x 2 \x-y\ +x 2 

\x-y\>- (|xi - 2/1 1 + F2 + 2/2 1 ) > - > - • 

Finally, we obtain the following estimate for L\. 



(2.9) 



\Li{x,y)\ < 



\x - y\(\x - 2/| + x 2 ) 



In order to estimate u 2 (t,x) we first estimate the contribution of vorticity-bearing par- 
ticles far from x: 



(2.10) 



'{|v-*l>i} 

From (|2.9j) and (|2.10[) we deduce that 



\L 1 {x,y)\u{t,y)dy<A / u(t,y)dy < 4\\u \\ L i. 



:2.11) 



\u2(t,x)\ < Cx 2 + Cx 2 



1 



{i«c-»i<i} \ x -y\\\ x -y\ + x 2. 



■u(t,y)dy 



for some constant C > 0. Changing to polar coordinates and estimating u(t, •) by the 
L°°-norm of Uq, it follows that 



f 1 1 

\u 2 (t,x)\ < Cx 2 + Cx 2 / — -r dr = Cx 2 + Cx 2 (log(l + x 2 ) - log 

Jo r(r + x 2 ) 



x 2 ), 



for some constant C. Relation (12.81) now follows. 



□ 



Remark 2.3. We call attention to the fact that this result holds without any condition on 
the sign of vorticity. 
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2.4. Vertical confinement. Another piece of information on confinement of vorticity for 
half-plane flows stems from the conservation of the second component of the center of 
vorticity. Such a result was obtained by Iftimie in [O] (see Theorem 3 and also Remark 3 



of (TT[). The resulting estimate will be used in what follows, so that we include its precise 



statement here for the sake of completeness. 

Theorem 2.3 ([0). If the initial vorticity ujq belongs to L P C (H), p > 2, then there exists 
a constant D > 0, depending solely on ujq and p, such that 

suppcu(t, •) C {x G H ; x 2 < D(tlogt)*} 
for all t > 2, where uj(t, ■) is a weak solution of having initial vorticity u . 



The complete proof can be found in 1 11 ] . 



3. Asymptotic behavior of nonnegative vorticity in the half-plane 

We now turn to our main concern in this paper, the rigorous study of the asymptotic 
behavior of flows with nonnegative vorticity in the half-plane. We divide this section in 
three subsections. In the first one we introduce the self-similar rescaling of the flow which 
encodes the scattering information we wish to study, we write an evolution equation for 
the rescaled vorticity and we interpret the vortex confinement information obtained in the 
previous section in terms of the new scaling. The second subsection is the technical heart of 
this article, where we study the behavior of the nonlinearity in the equations with respect 
to the self-similar scaling. Finally, in the third subsection we use the information obtained 
to prove our main result. 

3.1. Rescaled vorticity and asymptotic densities. One key feature of vortex dynam- 
ics in a half-plane is nonlinear wave propagation. In order to examine wave propagation it 
is natural to focus on a self-similar rescaling of physical space, as has been performed by 
Chen and Frid in the context of systems of conservation laws, see ||. Let us fix, through- 
out this section, a nonnegative function ujq £ p > 2, and uj = u(t, •), u = u(t, •), 
solutions of Ql.l|) with initial vorticity ujq. Set 

(3.1) uj(t,y) = t 2 uj(t,ty) and u(t , y) = tu(t , ty) , 

the rescaled vorticity and velocity, respectively. The scaling above respects the elliptic 
system relating velocity and vorticity so that we still have 

div u = 
curlw = uj. 

It is immediate that v,2(t, x±,0) = and therefore we can recover u from uj by means of 
the Biot-Savart law for the half-plane: 



(3.2) u(t,x)= I K(x,y)uj(t,y)dy, 

with K defined in (O). 
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Let M = \\u\\l°°cr + xw)- Then the confinement estimates for vorticity in the half-plane, 



in particular Theorems 2.1 and 2.3 and the fact that the vorticity u is transported by the 



velocity u, imply that there exists a constant C > such that: 

support, •) C [-C(t]ogt)*,C + Mt] x [0,(7(tlogt)5] for all t> 2, 

where Co = sup{xi ; x G suppc^o}- This in turn implies a sharp asymptotic localization 
on supp uj{t, •) = j support, •), namely: 



(3.3) suppS(V)C [-C( l ^,^ + M} x [0,C(^) ±3 



-logtx il 

Next, from the vorticity equation one may derive a transport equation for the evolution 
of u(t, y), which takes the form: 



(3.4) d t u{t, y)-- div [yu(t, y)] + - div [u(t, y) u(t, y)] = 0. 
Using the scaling (|3.1| ) we find 

(3.5) || u(t, OHw = t 2 ^\\u{t, -)\\ Lq < t 2(1 -^\\uj \\ Lq \/q E [l,p\. 

Furthermore, the L 1 -norm of uo is conserved in time. We wish to treat u as a bounded 
L 1 -valued function of time, possessing nonnegative measures as weak-* limits for large 
time. The confinement estimate ( |3.3|) implies that any weak-* limit of u must have the 
structure \i ® ^0(^2), with the support of fi contained in the interval [0, M}. 

It is in the nature of the self-similar rescaling ( [3.1|) that much of the scattering behavior 
of the flow is encoded in the measure \x. This measure is the main subject of the remainder 
of this article, and, as such, deserves an appropriate name. 

Definition 3.1. Let jj, G BM([0,M]) be a nonnegative measure such that there exists a 
sequence of times — > 00 for which 

u(tk, ■) — 51 fJL ® Sq in the weak-* topology of bounded measures, as t^ — ► 00. 

Then we call fi an asymptotic velocity density associated to ujq. 

It can be readily checked that, if u(t,x) = uj (xi — at,x 2 ), then there exists a unique 
asymptotic velocity density /i, which is a Dirac delta at position (a, 0) with mass given by 
the integral of lvq. For a general flow an asymptotic velocity density encodes information 
on typical velocities with which different portions of vorticity are traveling. 

3.2. The key estimate. Our purpose in this article is to understand the structure of the 
asymptotic velocity densities. To do so we make use of the evolution equation ( |3.4|) for co 
and we examine the behavior for large time of each of its terms. The main difficulty in 
doing so is understanding the behavior of the nonlinear term div (u u), which is our goal 
in this subsection. 
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We begin with two general measure-theoretical lemmas which will be needed in what 
follows. These are standard exercises in real analysis and we include the proofs only for 
the sake of completeness. Recall that a measure is called continuous if it attaches zero 
mass to points. 

Lemma 3.1. Let n be a finite and compactly supported nonnegative measure on R. Then 
fi is the sum of a nonnegative continuous measure v and a countable sum of positive Dirac 
measures (the discrete part of fi). Moreover, for every e > there exists 5 > such that, 
if I is an interval of length less than 5, then v(I) < e. 

Proof. Let A = {x ; fi({x}) ^ 0}. Then A is countable; indeed, A = [j n A n and each 
A n = {x ; fi({x}) > l/n} must be finite because \i is finite. Hence we may write A = 
{xi, X2, • • • } and rrij = fi({xj}). Of course, v — /i — £V rrij5 x . is a continuous, nonnegative 
measure. 

Let J be a compact interval containing the support of ji. For each x e J, it follows 
that, since v({x}) = 0, there exists S x > such that u(\x — 5 x ,x + 5 X ]) < e/2. Let 
Jx — [ x — &xi x + S x ]. Then J C {J x€j J x so that, using the fact that J is compact, we can 
extract a finite subcover, J C J Xl U • • • U J Xn . The interval J is now divided in a finite 
number of disjoint intervals (not necessarily J Xl , J X2 , . . . , J Xn ), each having //-measure less 
than e/2. It is then sufficient to choose 5 equal to one half of the minimum length of these 
intervals. For that choice of 5, it is clear that an interval of length 5 cannot intersect more 
than two of the disjoint intervals constructed above so that its //-measure will be less than 
e. □ 

Lemma 3.2. Let ^ n be a sequence of nonnegative Radon measures on M, converging weakly 
to some measure 7, and having the supports uniformly bounded in the vertical direction. 
Then, for every compact interval [a, b] one has that 

limsup7 n ([a, b] x R + ) < j([a, b] x R + ). 

n^oo 

Proof. Fix e > 0. Since 

7([a,6] x R+) = lim7([a - 5,b + 5] x R+J, 

8— >0 

there exists 5 > such that 

7([a - 5, b + S] x R + ) < 7([a, b] x R + ) + e. 

Let (f be a continuous function supported in (a — 5, b + 5) and such that < tp < 1 
and v|r ob ]— 1- According to the hypothesis, we have that (7 n (y), <p(yi))— > (7(2/), v(yi))> 
so there exists iV such that 

( 7n (y), ^(yi)>< (7(2/), ^(yi)>+^ Vn > N. 
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From the hypothesis on the test function cp it follows that, for all n > N, 

7«([a,6] x R+) < {in{y),^{yi))< {l(v)Mvxi)+ e 

< -/{[a - 5,b + 5] x R + ) + e < -/{[a, b] x R + ) + 2e. 

We deduce that 

limsup7 n ([a,6] x R + ) < 7([a,6] x K + ) + 2e. 

n^oo 

The desired conclusion follows by letting e — > 0. □ 

Let us now return to the study of the asymptotic behavior of vorticity. Let ujq > be a 
fixed function in L£(H), for some p > 2, and let it, u; be solutions of (jl.lp , with w, cD defined 
in ( [3.11) . Let /i be an asymptotic velocity density associated to Uq. Then /i is a nonnegative 
measure in BM([0, M]), with M = H^H^oo, and by Lemma |3TT|, fi can be written as 



(3.6) [L = y + y^ffm^, 

i=l 

where z/ is the continuous part of /i and ctj G [0,M]. As cj > it follows that m 8 > 
and, as is a bounded measure, m « < 00 • Furthermore we can assume without loss 
of generality that 7^ aij in the decomposition (|3~^ ). 

Let {tfc} be a sequence of times approaching infinity such that 

uj(t k , •) (1 ® ^0(0:2), 

as — > 00, weak-* in 5M(H). The following proposition is what we refer to as the key 
estimate in the title of this subsection. 

Proposition 3.1. Let ip e C (R). Then there exists a constant D > 0, depending only on 
p, such that the following estimate holds: 



(3.7) limsup 

fc— >oo 



f ^ yi )^iyl^(t k ,y)dy <D\\u \\if2 m T^\^i)\- 



Remark 3.1. It will be clear from the proof that the constant D can be chosen as D = 
Di^tt^ 1 , where Di p is the constant of Lemma f2.1| . 

Before giving the proof of Proposition |3.1|, let us motivate the statement with the follow- 



ing example. Consider a steady vortex pair with vorticity given by u(t, x) = ujq{x\ — at, x 2 ) 
and velocity u(t,x) = Uq(xi — crt, x 2 )- Then it is easy to see that the rescaled nonlinear 
term ^ uj converges to amb a ® 5q where m = J u dx. Based on this example, one would 
expect the right-hand side of ( p.7|) to be J2i a i m i\' l l } ( a i)\ instead. On the other hand, for 
the steady vortex pair, it can be easily checked that 

a = 1 / (M )!u;odx < |K||i/*>- 
J uj dx J 
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Using Lemma [2.1| we infer that 



P'/2 l-p'/2 



\a\ < D\\uJ Q \\ L p m 

which then implies that, as measures, the weak limit of ^- u is less than D\\ujQ\\ p L l' 2 'm 2 ~' p '/ 2 5 a ® 
5q. Hence, in light of this example we see that estimate ( |3.7|) is weaker than what might 
be expected, but nevertheless it is consistent with the behavior of steady vortex pairs. 



Proof of Proposition \3. 4 Let us denote the integral we wish to estimate by B^, so that 



(3.8) 



-Lo(t k ,y)dy. 



Fix e > throughout. Since Y^=x m « < 00 there exists iV = N(e) such that 



E 

i>N 



e 

m,i < -. 
4 



(3.9) 



Additionally, it is easy to find 5 = 5(e) > such that, if / is an interval, |/| < 5, then 

e 



!/(/)< 



l,...,iV, 



by using Lemma [3.1| , and also 

(3.10) fi([ai-26,ai + 26\)<rm0- + e), 

(3.11) [on - 5, oa + 8] n [aj -5,a j + 5}=dl 1 i ^ j e {1, . . . , N}, 

(3.12) \ij)(y 1 )-^(a i )\<e V ^ e [a, - 5, ^ + 5), i = 1, . . . , AT. 

In view of Lemma |3.2| and relation ( |3.10 ), there exists K such that, if > K then 



(3.13) 



/ 



uj(t k ,y)dy < m;(l +e) 



Vi = 1, . . . , N. 



[ai-2<5,ai+2<5]xIR + 



N 



Consider now an interval I C R \ |J (at* 



relation ( PD 



|, «j + |) of length at most 5. According to 



u(I)< 



On the other hand // — u, the discrete part of /i, restricted to / avoids the Diracs at 
cti, . . . ,ajsr so that 

Therefore 

(3.14) /i(J) < |. 
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N 

Given a compact interval J C R \ |J (ol^ — |, aij + |) of length at most 5 we can use 

i=l 



( |3.14| ) and Lemma |372| together with the fact that ui(t k , •) — ^ jj, (g><5 to find iT large enough 

w(t fc ,y) dy < 



so that, in addition to ( |3.13|) , we have 

e 



9 

for any k > K Q . We wish to show that this K can be chosen independently of 3 , but we 
shall have to pay a price, namely the estimate above will hold with e on the right-hand-side, 
instead of e/2. 

Let J be a compact interval such that J x R + contains the support of u(t, •) for all t. 

N 

We write the set J \ \J (cti — |, ctj + |) as a finite disjoint union of intervals Ij, each of 

i=l 

which we subdivide into intervals of length exactly 5, together with an interval of size at 

N 

most 5, this being the right-most subinterval of Ij. This way the set J\ \J (on — | , on + f ) 

i=i 

can be written as the union of intervals J\ , . . . , Ji of length precisely 5 plus some remaining 
intervals Ji + i, . . . , Jl of length strictly less than 5. According to ( |3.14| ), we have that 

//(Ji) < | Vz = 1, . . . , L. 



Next we apply Lemma 3^2 and use the fact that u(t, •) — ^ jj <8> 5, to obtain i-To such that 



( p,13| ) is satisfied together with: 

(3.15) J u(t k ,y)dy< S - Vz = l,...,L, k > K . 

N 

Let / be a subinterval of R\ IJ (at* — |, + |) of length less than 5. It is easy to see that 

i=l 

I can intersect at most two of the intervals Jj as otherwise, by construction, this would 
imply it had to contain an interval of length precisely 5. According to ([3.15|) we deduce 
that fj M u(tk,y) dy < e for all k > K Q . We have just shown that, if I is an interval of 

N 

length at most 5, I C R \ (J (aij — f , ctj + f ) then 

(3.16) J u(t k ,y)dy<e, Mk > K . 

IxR + 

Let k > K and set 

E t = [ai - 5, at + 5] x R + , Fi = [oti - 25, on + 28] x R + , £7 = E x U • • • U E N . 
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According to ( |3.11| ), the sets Ei,...,En are disjoint, so we can write B k , defined in 
TSp , as: 

i=l J E i k J Ec _ fc , 



^ „ V 

Bh.o. 



B k i 

We will estimate separately Bki and B k 2- Note that both estimates rely in an essential 
way on the Biot-Savart law and the fact that the kernel can be estimated by \x — y\~ x (see 
( p.4| )). In the remainder of this proof we will denote by C a constant which is independent 
of e and t. 



Estimate of B k \. Using the Biot-Savart law ( |3.2| ) and relation ( |2.4| ), one can bound 
Bki as follows: 

\B»\<tf[ l f ( ' ft)l | g( V , 'W »)d*d» 

I I n\x-y\ t k 



1-1 zee 

N 



r/. // —\ :uj{t k ,x)uj(t k ,y) dxdy 

'kT-f J J n\x-y\ 



|x— y|><5 



+ 7-^ // f . u(t k ,x)u(t k ,y)dxdy 
tk~ J J n\x-yi 



\x— y\<S 



j=l j=l 

SD— 3/| <0 

rreH,j/e£; 



^JT + ]C // I &(t k , x)uj(tk, y) dx dy. 



% \x-y\<& 



According to (|3.12|) , for y e Ei we have that |^(yi) — if>(oii)\ < e - We therefore deduce 
that 

\x-y\<6 
x£U,y£Ei 
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Applying Lemma |2.1| yields 



u(t k ,x)u{t k ,y)dxdy < [ ( I ( £_^!L d x )u(t k , y) dy 



\x-y\ Je^ J \x-y\ 

\x-y\<8 [j/i xR+ 



< D 1>p J ( J u(t k ,x)dx S j 2 || u\\l P u(t k ,y)dy. 



[y 1 -S,y 1 +S]xR+ 

Now, if y e Ei then [yi - 8, y 1 + 5} C [a { - 25, a, + 25], so that [yi - 5, y 1 + <5] x R + c F { . 
Hence 

j-^—^(t k ,x)u(t k ,y)dxdy<D hp ^J u{t k) x)dx s j * t k \\uJo\\% (^J u(t k ,y)dy 

\x-y\<8 
xeM,y£Ei 

<t k D ltP [ mi (l + e)] 2 ~ P ~\\uj \\%, 
where we have used (|3.5| ) and (|3.13|) . We conclude that 

(3.17) [-Sfcil <— + C 1 J2(m^)\+e)[m l (l + e)] 2 -^, 

° tk i=i 

eL 

with d = Di^HwollipTT \ 

Estimate of l?^. We estimate directly, similarly to what was done with B k \: 

\B k2 \<£- + II J^ (Z/l)l r u(t k ,x)u(t k ,y)dxdy 
ot k J J ■Kt k \x-y\ 

\x-y\<5/Z 
xm,y£E c 



- TT + ^^] L °° // ] — - — r a;(ti fc ,x)a;(t fe ,y)dsdy. 

\x—y\<8/3 
xeU,y£E c 

Lemma [2.1| implies in the same way that 

l-Bfcal < t— + — r^lMU 00 / ( / u(t k ,x)dx) 2 \\uj\\I p u(t k ,y)dy 

ot k irt k J E c\ J / 



CD — 

<— H ^||^|| L oo||u; || L 2 p ||^ || L i sup ( / u(t k ,x)dx 

0t k 7T yeE < 



1-2- 



[2/1-1.2/1+1] 
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JV 

For y e E c , the interval [y% — |, y% + 1] is of length less than S and included in 1R\ |J (ctj • 

1=1 



§, OLi + |). We deduce from ( |3TT6l ) that 

u3(tfe, x) dx < e. 



/ 



[yi-f,ltt+f]xR+ 

which implies that 

(3.18) |5 fe2 | < ^ + C 2 e 1 -^ 



with C 2 = C , i||V ; IU° o ll^o|U 1 - 



Collecting the estimates for Bki and B^ (relations (|3.17|) and (|3.18|) ) yields the following 
bound for B^. 

n N ' 

(3.19) \B k \ < — + C 1 {||^|| L oo|| W o||L 1 ^^ + ^(|^(«,)|+^)h^( 1 + £ )] 2_ ^}- 

Take the lim sup as k — > oo above to obtain: 

oo , 

limsup|£ fe | < d{ II^IUooll^olUi^ 1- ^" +53(1^)1 +e)[m t {l + e)] 2 ~ P ~\. 

1 = 1 

Next, send e — > in order to reach the desired conclusion. □ 
3.3. Large time asymptotics. We will now make use of the equation for uj given in 



relation (p.4[) together with Proposition 3.1 to deduce an inequality for the limit measure 



\x 1 given by ( ^.23[ ). Surprisingly, this estimate alone will be sufficient to deduce the main 



result of this paper, Theorem |3.1| . Let us begin with an outline of the proof of ( |3.23j ). 



One begins with the equation for the evolution for uj ( |3.4| ), taking the product with a 
fixed test function and integrating in space. The resulting equation has three terms. The 
first one, when integrated from to t, is uniformly bounded in t. Now, if div[y uj(t, y)] is 
weakly convergent as t — ► oo, then the integral in time of the second term will, in principle, 
diverge like \ogt as t — > oo. As for the third term, it is not difficult to see that it is 0(l/t). 
The dominant part of the third term must balance the logarithmic blow-up in time of the 



second term. The aim of Proposition |3.1| is precisely to estimate this dominant part of the 
third term. 

We will begin with a lemma, relating asymptotics on the linear part of the evolution 
equation for uj (|3.4|) to the nonlinear part. To this end fix ip G C°(R) and define the 
quantities 

(3.20) A[t;iP}= [ iJ>(yi)yMt,y)dy, and 
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(3.21) B[t-M = [ ^ yi )^ylz(t,y)dy. 

Jw 1 

Note that, as the support of uj is contained in a compact set independent of t, it will not 
matter whether the support of ip is compact. 

Lemma 3.3. The following estimate holds: 

(3.22) limsup(£[t; tp) - A[t; ip)) > 0. 

t— >oo 

Proof. Let <p G C 1 (1R) be a primitive of ip so that cp' = ip. Define 



/(*) = / f(yi) v) d v> 

Jw 

a bounded function, since u(t, •) is bounded in L 1 . Differentiating / with respect to t and 
using the equation ( |3.4j) for u we get, after integration by parts, 

f(t)=J <p(yi)d t ti(t,y) dy=jj (p(yi)eiy(yu) dy - ^ j p( yi ) div(uuj) dy 

1 f XL 1 C 1 1 

= j / V'Cl/Oywdy- - ip(y 1 )y 1 udy = -B[t;ijj\--A[t;'ip]. 
Integrating from t to t 2 we obtain: 

Let L = limsup s _ +oc (-B[s; ^] — ^[ s lV'])- Then, for any e > 0, there exists M > such 
that, if s > M then S[s; ip] ~ + In particular, if t > M above then 

f(t 2 )-f(t)<(L + e)\ogt, 

so that 

i^oo log t 

The result follows by taking e — > 0. □ 

Let us now impose a major hypothesis on the flow, namely that there exists a unique 
asymptotic velocity density, so that 

uj(t, •) \i ® <5(x 2 ), 



as t — > oo. We use Lemma |3.1| to write 

oo 



Then, for any ^ G C°(R), it follows that 

A[t;i/j] -> (yifi,ip(yi)), 
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as t — > oo. Next use Proposition [3J] to deduce that 

, oo , 
E- \ -» 2-2- 

limsup \B(t; ip}\ < D\\u \\l P m i 2 I^O**)!- 



i=i 

We therefore deduce from Lemma RT31 that: 



i=i 

Exchanging ip for —ip yields: 

(3.23) |<2/i^V>(2/i)>| <D\\uo\\%J2 m ^^ a ^- 



2-8- 

rr 

8=1 



The relevant fact is that the exponent 2 — ^ > 1. 

Let <5p denote the Dirac delta measure at position P. We are now ready to re-state our 



main result, giving a more precise formulation of Theorem 1.1 



Theorem 3.1. Suppose that the nonnegative initial vorticity uq G L p c {M) is such that there 
exists a unique asymptotic velocity density fi associated to ujq. Then [i must be of the form: 

oo 
i=l 

where: 

(a) cti 7^ acj if ' i ^ j and ai — > as i — > oo; 

(b) the masses rrii are nonnegative and verify YliLi m i = ll^olU 1 / 

(c) for alii, a { e [0,M], where M = ||n|| L oc( [0>oo)xH) ; 

(d) there exists a constant D > 0, depending solely on p, such that, for all i with rrii ^ 
we have 

E- 1—2- 

ati < D\\u \\l P m i 2 . 

Furthermore, there exists i^ such that a io ^ and m io ^ 0. 

We will need two lemmas before we give the proof of Theorem |3.1| . 

Lemma 3.4. Let mi5 ai be a Dirac from the discrete part of fi. The following inequality 
holds true: 

cti < D\\u \\l P m i 2 . 

Proof. Eventually changing the order in the summation of the Diracs, we can assume that 
i = 1. Furthermore, the conclusion is trivial if a± = 0, so we can assume that a\ > as 
well. Let e > be fixed. There exists S e (0, «i) such that the following inequality holds: 

fi([ai — S, ax + S\) <mi+e. 
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If aj G [cki — 5, a>i + 8], i > 2, then m\8 ai + mi5 ai < /i on [a% — 8, ai + 5], so we must 
have that rrii < e. 

Let ^ G C°(R) be a nonnegative function supported in (ai — 5, «i + 5) C R+ which 
attains its maximum at a.\. By fl3.23| ) and using the nonnegativity of \i and yiip(yi) we 



find 

P ■ > ; ■■ — — e, ■ > j 

miai^(ai) < (n,yitp(yi))< £>||a;o|| L 2 p [^("1)^1 2 + 2^ ^(a i )m i 2 ] . 



i=2 



We observed that if G («i — 5, a\ + 5), i > 2, then raj < e. If G" («i — 5, «i + 5) 
then ip(oii) = 0. In both cases 



r tl){a.,i)rn i 2 < tp^a^e 1 " 2 m,i < ipia^e 1 ^ 
We infer that 



r 2 — i V V — v n 

miaiV'Cai) < Z?||a>Q||ipV(ai) L m i 2 + £ ~ T / 

1=2 

that is 

/ / .00 

P r- 9—P 
I O ^ 



■m\a\ < D\\uj \\l P [m 1 2 + e 2 \ mi j. 



i=2 



Letting £->0we get that 



'! 2 



m\a.\ < D\\ujQ\\l P m 1 

which implies the desired result. □ 
Lemma 3.5. Suppose that [i has no discrete part in some interval (a, b) C M\ {0}. Then 

Proof. Let if> G C°(1R) with support in (a,b). According to the hypothesis, 

supp if) fl {ax, «2, . . . } = 
so that, for this choice of if), the right-hand side of ( |3.23| ) vanishes. Therefore ( |3.23| ) implies 

(v(yi),yiip(yi))= o 

that is 

which implies the desired conclusion by recalling that G" (a, b). □ 
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Proof of Theorem We begin by noting that Lemma pT4| implies that cti t ^°°> 0. Indeed, 

X^i m « < 00 implies that mi ^°°> 0. According to the conclusion of Lemma |3.4j this 

immediately implies that a, t ~ > °°> 0. 

Next, observe that Lemma [3.5| implies that the continuous part v vanishes. Indeed, 
supp v C [0, 00) since supp ji C [0, M]. If a > 0, as a is not an accumulation point of the 
set {ax, a 2 , ■ ■ ■ }, there exists S G (0, a) such that {ax, ct2, . . . } PI [(a — 5, a + 5) \ {a}] = 0. 
According to Lemma |3.5|, the measure /1 vanishes in (a — 5, a) and (a, a + 5), so the same 



is true for v. Since v is continuous we deduce that v must vanish in (a — 5, a + 5). We 
proved that v vanishes in the neighborhood of each point of (0,oo). This implies that v 
vanishes on (0, 00). Therefore, v vanishes on R \ {0} and is continuous. We conclude that 
v = 0. 



We have just proved that 



i=l 



and also assertion (|a|) of Theorem |57I . Assertion (|Bp follows from the positivity of /x (as 
limit of positive measures) and from the fact that the total mass of \i is H^ollz, 1 - Assertion 
(H) is a consequence of the support of \i being included in [0, M\ and (^) is proved in 
Lemma ^.4| . Finally, as previously noted, it was shown in |TT| that J xiu(t,x) dx > Ct 
for some positive constant C. This implies that J xiu(t, x)dx > C, which in turn yields 
J2i m i a i — (a*) C . This completes the proof of Theorem □ 

4. Extensions and Conclusions 

We begin this section with some comments regarding the results obtained here. 

(a) The only instance of use of the energy estimate in this work is the observation that, 
for any asymptotic velocity density /x we have (/x, Xi)> C > 0, which appears when 
proving the last part of Theorem |3.1| . The constant C depends on the kinetic energy 
of the initial data, as was derived in [IT]. It would be interesting to know whether 



kinetic energy partitions itself in a way that is consistent with the partitioning of 
vorticity, but we were not able to prove that, at least using only the hypothesis of 
uniqueness of the asymptotic velocity density. 

(b) We only used the hypothesis of uniqueness of the asymptotic velocity density when 
we derived ( [3. 23] ). The estimate on the behavior of the nonlinear term given in 
Proposition |3J] always holds, which raises the possibility of it being exploited further. 

(c) The hypothesis that the initial vorticity be p-integrable, with p > 2 is used to ensure 
that the velocity is globally bounded. In principle, with vorticity in L p , p < 2, we 
loose control over the loss of vorticity to infinity, and Lemma [2.1| is no longer true. In 
fact, we do not even know the correct scaling to analyze in this case. 

(d) We could have proved a result for initial vorticities without distinguished sign, which 
would have our main theorem here as a special case. The distinguished sign hypothesis 
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was used mainly because it allows us to place the hypothesis of uniqueness of the 
asymptotic limit on cu, instead of on \uj\, where the authors feel it would be much less 
plausible. 

We would like to add a remark on the choice of the scaling x = ty. If the scaling 
x\ = tyi in the horizontal direction is motivated by the fact that the first component 
of the center of vorticity behaves exactly like 0(t), the scaling x 2 = ty 2 is not justified 
because the second component of the center of vorticity is constant. Ideally we should 
not make any rescaling in the vertical direction but then we would have to assume that 
tu(txi, x 2 ) converges weakly which we found excessive because of the oscillations that 
may appear in the vertical direction. We could also consider an intermediate scaling of 
the form x 2 = f(t)y 2 where f(t) — > oo as t — > oo. This last problem is in fact equivalent 
to the one we consider in this paper. If / is such a function, then the weak limits of 
cijf(t,y) = tf(t)uj(t, tyi, f(t)y 2 ) are independent of /. Indeed, let Vf be the weak limit of 
ujf(t,y) as t — > oo and choose a test function h G C^°(H). Then 

( z//) h) = lim / u f (t, y)h(y) dy 

= lim / u{t,x)h(^,^-)dx 

— lim ( [ uj(t,x)h(^,0) dx + o(—^-^—) [ x 2 u(t, x) dx) 
*^°° v Jh t V }{t) J J m J 

= lim / uj{t,x)h{ — , 0) dx 

since we know that f m x 2 u(t,x) dx = est. and f(t) — > oo as t — > oo. The last term does 
not depend on / anymore. Here, we have made the choice f(t) = t only for the sake of 
simplicity. This means that we study the asymptotic behavior of solutions in horizontal 
direction but not in the vertical one. 

We would like to comment on a few problems that arise naturally from the work presented 
here. The first is to remove the hypothesis of uniqueness of the asymptotic velocity profile, 
perhaps with weaker conclusions. Also, we can try to extend this line of reasoning to other 
fluid dynamical situations with similar geometry, such as flow on an infinite flat channel, 
axisymmetric flow (smoke ring dynamics), and water wave problems. We may also ask the 
same questions with respect to full two-dimensional scattering, allowing for vortex pairs 
moving off to infinity in different directions. Finally, one might try to examine the issue 
of actually proving the uniqueness of asymptotic velocity densities in special cases, for 
example, for point vortex dynamics. We have actually looked at the case of three point 
vortices on the half-plane, so far without success. 

Appendix. Separation of two vortices above a flat wall 

Steady vortex pairs provide smooth examples of vorticities for which the corresponding 
asymptotic velocity densities consist of a single Dirac mass. We would like to give such an 
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example with at least two different Dirac masses in the asymptotic velocity density. As we 
pointed out in the introduction, the existence of multibump solutions in this situation is 
an interesting open problem, but we can offer a discrete example in order to illustrate this 
issue. In this section we will give a sufficient condition for linear separation of two vortices 
above a flat wall which will in turn give us an example of unique asymptotic velocity 
density concentrating at two distinct Dirac masses. 

Let z\ = (xi,yi) and z 2 = (£2,2/2) be two vortices above the wall {y = 0} of positive 
masses 777,1, resp. 777.2. For notational convenience we will assume that we start at time 
t = 1 instead of t = 0. Let L be defined by 

(A.l) L = 777,12/1 + 777 2 7/2, 

a quantity which is conserved by the motion of the vortices. We will prove the following 
proposition. 

Proposition A.l. Suppose there exists a positive constant M such that the following re- 
lations hold true: 

(A.2) x 2 (l) - rri(l) > M, 

(A.3) L > m 2 y 2 (l) + ^ 

and 

(A.4) - - ^ ^ - 2max(m 1 ,m 2 ) > 2nM. 

Then, the two vortices Z\ and z 2 linearly separate. More precisely, 
(A.5) x 2 {t) -xiit) > Mt 

for all times t > 1 . 

Remark A.l. Let mi, m 2 and L be some fixed arbitrary positive constants. Then we can 
always find £i(l), 2/1(1), £2(1), 2/2(1) and M such that relations ( |A.1| ), ( |A.2|) , ( |A.3|) and 
( |A.4| ) are satisfied. Indeed, we first choose Xi(l) and x 2 (l) such that ( |A.2| ) holds. We next 
note that ( |A.3| ) and ( |A.4| ) are satisfied for large enough M and small enough 2/2(1)- For 



example, if 2/2(1) = 0, then ( |A.4| ) has a left-hand side of order M 3 so it is verified for M 
large enough; and since, for that choice of M, it is satisfied for 2/2(1) = 0, it will be satisfied 
for small enough 2/2(1), too. Once 2/2(1) and M are chosen, it remains to choose 2/1 (1) such 
that (|A.1|) is satisfied for t = 1. 



Proof of Proposition PT7| . It is sufficient to prove that, as long as ( |A.5| ) holds, then 



(A.6) (x 2 - zi)'(t) > M. 
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Indeed, the result then follows by a contradiction argument: if T is the first time when 
^2 GO — x i(T) = MT, then necessarily T > 1 and 

MT = (x 2 - x x )(T) = x 2 (l) - xi(l) + J (x 2 - an)' > M + M(T — 1) = MT 
which is a contradiction. 

We will therefore assume in the following that ( |A.5|) holds and try to prove (|A.6| ). 

It follows from the method of images that the motion of these vortices can be computed 
from the full plane flow due to these two vortices together with their images: 

z 3 = zi = (xi, ) and z 4 = z 2 = (x 2 , -y 2 ) 

with masses m 3 = —mi, resp. = —m 2 . Therefore, the equations of motion are given 
by: 

„ / Oi - Z2) 1 - . {zi - z 3 ) x {zi - z 4 ) X 

2nz 1 = w m 2 + -, nrm 3 + pr m 4 ; 

\z 1 -z 2 \ 2 \z 1 -z 3 \ 2 \z 1 -z 4 \ 2 

i.e., 

2Trz' 1 = 2Tr(x' l7 y' 1 ) 
(A.7) ( m i n\ m 2 / m 2 

2?/! |^i -^r 1^1 -22f 

Interchanging the indexes 1 and 2 we also get 

2ttz 2 = 27r(x 2 ,y 2 ) 

(A.8) ,m 2 , m 1 , mi , 

= (tt - ' °) + i n(yi - 2/2, 2:2 - an) + , =-r^U/i + 2/2, an - x 2 ). 

y 2y 2 \zx-z 2 \ l \z\-z 2 \ l 



Let us now estimate y 2 . From relation (|A.8| ) it follows that 

1 1 \ mi(x 2 - x 1 )Ay 1 y 2 



27ry 2 = mi(i 3 - an) f, ^ - , 

\\z-\ — zo\ z \z-i — Zor/ 



Z\ — Z 2 \* \Z\ — z 2 \* / \z\ — z 2 \ 2 \zi — z 2 \ 2 



In view of ( |A.lj ), we can bound miyi < L and y 2 < L/m 2 so that, using also relation 



( A -9) Va < — TZ ^T, < 



nm 2 \xi — x 2 \ 3 nm 2 M 3 t 3 



We deduce that 



Mt) - y 2 (l)\ = I y' 2 <-^r 3 f 4d S = -^(l-4)< 



which implies that 

(A.10) j, a (t) < t/ 2 (l) + 



7rm 2 M 3 J 1 s 3 nm 2 M 3 t 2 nm 2 M 3 
L 2 



irm 2 M 3 
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Next, from (|A.7|) , (|A.8|) and (|A.1|) we have that 



(x 2 - x x )' 



> 



1 


~m 2 


m 1 | (mi + m 2 )(2/i -2/2) 
2yi ki-^| 2 


(mi 


- m 2 ){yi+y 2 ) 


2^ 


-2^ 




\zi -z 2 \ 2 


1 


"m 2 


m 2 (mi + m 2 ) 


mi 


- m 2 1" 




2^ 


-2m 


2(L-m 2 y 2 ) \z\-zr\ 




-^2 - 





Both |zi — z 2 1 and — ~z 2 \ are bounded from below by \x± — x 2 \ > Mt > M. Furthermore, 
the first two terms of the right-hand side of the last relation are decreasing with respect 
to y 2 . We therefore deduce from ( |A.10| ) that 



(x 2 - Xi)' > 



1 

2^ 
1 

~ 2^ 
> M, 



m 2 



(mi + m 2 ) 



mi — m 2 



m 2 



mt 



2(|/2(1) 



L 2 



) 2(L-m 2 2/ 2 (i; 



L 2 ' 

7rM 3 , 



M M 

2max(mi, m 2 )" 
M . 



where we have used (|A.4j ). This completes the proof. 



□ 



Remark A. 2. The conclusion that x 2 (t) —X\(t) > Mt, for some M > 0, always implies the 
existence of a unique asymptotic velocity density which concentrates on a pair of Dirac 
masses. In order to see this, first note that, from (|A.9| ), we have that \y' 2 \ = 0(l/t 3 ), 
which implies that y 2 {t) converges as t — > 00 and similarly for y\. From the conservation 
of energy we have that 



2m\m 2 log 



\z\ - z 2 \ 



- m\ log(2yi) - m\ log(2j/ 2 ) 



is constant in time. Since x 2 {t) — X\{t) > Mt we also know that Z2 ^ 



\z\ - z 2 \ 

1 as t - 

We deduce that lim y 2 (t) 7^ and lim y\{t) 7^ 0. Now, from relations ( |A.7| ) and (|A~i 
immediately obtain that both x[ and x' 2 converge to a finite limit given by 

m l 1 T / /A m 2 



OO. 

we 



lim x[(t) 

t— >oo 



47r lim yi(t) 

t— >oo 



and 



a 2 



lim x' 2 (t) 



Att lim y 2 (t) 

t—>oo 



Observe next that lim = lim x'^t) = a± and similarly for 5^SL. Finally, let us remark 



^2 L^j — lim ^"/"^ — T- o n ciTnilarlif fjr\T" "^ 2 
t— »oo * t^OO 

that the rescaled vorticity is given in this case by mi5 Zl / t + m 2 5 Z2 / t so that it clearly 
converges weakly to (mi5 ai + m 2 5 a . 2 ) <g> 5 . Moreover, x 2 (t) — Xi(t) > Mt implies that 
a 2 - ai > M > 0. 
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